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ABSTRACT 


Procedures  are  presented  whereby  natural  vibration 
and  transient  response  analyses  of  partially  submerged, 
supercavitating  marine  propellers  may  be  calculated  by 
means  of  a  standard  NASTRAN  structural  analysis  program. 

The  use  of  these  procedures  was  demonstrated  by  analyses 
of  an  existing  propeller,  for  which  some  experimental  data 
were  available.  Results  were  comparable  only  in  a  limited 
way,  but  peak  calculated  dynamic  stresses  were  of  the  same 
order  of  magnitude  as  those  observed  experimentally,  and 
were  apparently  more  than  twice  the  magnitude  of  equivalent 
static  stresses . 


ADMINISTRATIVE  INKORMATION 

Funds  for  this  investigation  were  obtained  from  the  Naval  Material 
Command  (NAVMAT  0333)  in  support  of  an  ongoing  Ship  Performance  and  Hydro¬ 
mechanics  Exploratory  Development  Program  (Program  Element  62543N,  Task 
Area  ZP-A3-421-0U1)  assigned  to  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC).  The  High  Speed  Propulsor  Task  of  this 
Exploratory  Development  Program  (DTNSRDC  Work  Unit  No.  1500-200)  provided 
funds  to  complete  work  presented  in  this  report. 

INTRODUCTION 

The  utilization  of  high  speed  marine  vehicles  by  the  United  States 
Navy  has  created  a  need  for  an  improved  structural  analysis  and  design 
capability  for  high  speed  propulsors.  These  are  generally  supercavitating 
propellers,  which  may  also  be  partially  submerged.  In  their  hydrodynamic 
performance  and  structural  integrity,  tliese  propellers,  because  of  their 
different  geometries  and  modes  of  operation,  represent  substantial  depar¬ 
ture  from  conventional  propellers,  for  which  a  large  body  of  experience 
exists.  Dynamic  aspects  of  structural  response,  in  particular,  are  ex¬ 
pected  to  be  more  significant  tor  these  propellers.  These  aspects  need  to 
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be  evaluated  so  that  attempts  to  attain  optimum  hydrodynamic  designs  will 
not  be  limited  unnecessarily  by  structural  considerations.  An  effort  was 
therefore  Instituted  to  develop  methods  of  predicting  dynamic  stress  and 
deflection  behavior  of  high  speed  propellers  under  realistic  time-dependent 
loading  conditions. 

A  preliminary  Investigation  Indicated  that  meaningful  results  were 
not  likely  to  be  obtained  other  than  by  computer  studies.  Various  computer 
programs  have  been  used  to  calculate  the  full  vibration  modes  and  frequen¬ 
cies  of  completely  submerged  propellers,  but  until  recently  the  effect  of 
Immersion  in  water  has  been  accounted  for  by  assuming  that  a  more  or  less 
arbitrary  mass  of  fluid  Is  attached  to  the  propeller  blade.  This  "added 
effective  mass"  approach,  which  Incorporates  the  assumption  that  the  fluid 
is  Incompressible,  is  considered  to  account  adequately  for  fluid  structure 
Interaction  in  most  vibration  and  transient  response  problems,  provided 
valid  added  fluid  mass  terms  can  be  obtained.  A  few  Investigators  have 
recently  published  descriptions  of  analyses  of  dynamic  behavior  of  sub¬ 
merged  structures  In  which  effects  of  submergence  in  water  are  computed  on 
a  rational  basis.  One  such  Investigation^*  yielded  a  procedure  by  which 
the  existing  NASTRAN  finite  element  structural  analysis  program  could  be 
applied  to  the  submerged  structure  problem,  even  without  making  the  incom¬ 
pressibility  assumption.  That  work  forms  the  basis  of  the  analyses  de¬ 
scribed  here.  The  approach  was  selected  because  of  the  availability  of 
NASTRAN  at  DTNSRDC  and  because  of  Its  flexibility  of  application. 


*A  complete  listing  of  references  is  given  on  page  57. 
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Further  modifications  were  made  to  this  procedure  to  adapt  it  to  the 
general  problem  of  predicting  free  vibration  frequencies  and  transient 
response  of  supercavitating ,  partially  submerged  propellers.  The  resulting 
approach  was  applied  to  determine  natural  frequencies  of  two  propellers, 
and  the  dynamic  stress  response  of  one  propeller  for  which  some  experimen¬ 
tal  data  were  available.  Only  a  limited  comparison  of  results  was  possible; 
however,  the  calculated  frequencies  and  dynamic  stresses  appeared  to  be 
reasonable.  The  dynamic  analysis  and  its  application  via  NASTRAN  are 
described  in  this  report. 


DESCRIPTION  OF  ANALYSIS 

In  the  analysis  to  be  used  here,  the  fluid  surrounding  the  submerged 
structure,  in  this  case  an  individual  propeller  blade,  is  represented  by  a 
finite  element  idealization.  More  precisely,  a  quantity  of  fluid  suffi¬ 
ciently  large  that  its  effect  on  the  structure  should  nearly  be  the  same 
as  that  of  an  infinite  quantity  is  represented.  The  theory  governing  this 
approach  is  described  by  Zienkiewicz A  concise  presentation,  with  a 
delineation  of  the  strategy  of  application  by  means  of  existing  NASTRAN 
capabilities  is  given  by  Everstine  et  al.^  The  theory  and  application 
procedure  will  be  briefly  summarized  here. 

Conventional  analytical  techniques  lead  to  the  following  matrix  equa¬ 
tion,  which  describes  the  dynamic  behavior  of  the  propeller  blade  when 
represented  in  finite  element  form: 

[mHu!  +  [kHuI  =  IF}  +  [Llipf  (1) 
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where  {u}  is  the  set  of  unknown  grid  point  displacements  and  is  a  set 

of  applied  time-dependent  forces.  The  column  matrix  fp}  is  the  set  of 
unknown  pressures  on  the  fluid-structure  interface  grid  points.  The  last 
term,  [L]{p}  >  represents  the  fluid  loading  on  the  structure. 

The  set  of  unknown  pressures  at  ail  the  grid  points  of  the  fluid  fi¬ 
nite  element  idealization  is  represented  in  the  following  matrix  equation: 

[  Q  1  \  p  ^  +  [  H  ]  i  p  ■  =  -  [  S  1  \  u  ■  (  2  ) 

The  last  term  in  Equation  (2)  represents  the  effect  of  accelera t ions  of  the 
structural  grid  points  at  the  fluid  structure  interface  on  Che  pressures 
in  the  fluid.  The  bar  placed  over  iu  •  indicates  that  only  normal  interface 
deflections  are  involved  In  that  terra.  The  fluid  raay  be  considered  an 
acoustic  medium,  in  which  case  Equation  (2)  provides  a  valid  description 
of  its  dynamic  behavior  if  it  actually  represents  the  wave  equation 


with  appropriate  boundary  conditions.  In  Equation  (3),  c  represents  the 
speed  of  sound  in  the  fluid.  It  was  shown  previously^  that  the  implemen¬ 
tation  of  the  procedure  can  be  accomplished  by  means  of  a  standard  struc¬ 
tural  program.  The  left  side  of  Equation  (2)  is  treated  as  being  analogous 
to  the  left  side  of  Equation  (1),  i.e.,  it  is  treated  as  a  dynamics  problem 
with  unknown  variables  Ip'  .  Only  one  translational  degree  of  freedom  at 
each  grid  point  of  the  fluid  is  left  unconstrained,  liowever,  and  that  one 
remaining  component  is  interpreted  as  pressure.  In  addition,  all  degrees 
of  freedom  may  be  constrained  on  the  outer  fluid  boundary  (ail  boundary 
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points  not  in  contact  with  the  moveable  structure).  Finally  for  the 
three-dimensional  case,  the  following  material  property  matrix,  repre¬ 
senting  Hooke's  law,  is  used: 
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With  the  NASTRAN  program,  it  is  not  usually  necessary  to  input  the  matrix 

3 

(4).  As  pointed  out  by  Everstine,  the  same  information  can  be  input  by 
specifying  an  isotropic  material  with  the  following  values: 

Shear  modulus  (J  =  1  .U 

Poissons  ratio  \'  is  not  input  (so  ttie  program  will 
compute  a  value). 

Young's  modulus  E  =  aO;  »  •'  1.0 

Here,  a+  l.U  should  be  numerically  indistinguishable  from  >.  On  most 
computers,  a  =  10^*^  will  suffice.  Equations  (1)  and  (2)  may  now  be  com¬ 
bined  into  one  set: 
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A  dynamic  problem  governed  by  Equation  (5)  can  be  run  (tlie  free  vibra¬ 
tion  mode  analysis  results  from  letting  F  be  null)  by  providing  finite  ele¬ 
ment  data  for  both  structure  and  fluid,  in  wiiich  case  all  matrices  except 


S  and  L  in  Equafo 


(S)  Will  be  set  up  automatically.  Note  that  grid 


points  for  structure  and  fluid,  including  those  on  the  interface  which  have 
the  same  coordinates,  are  numbered  separately.  Thus,  the  only  coupling  be¬ 
tween  fluid  and  structure  is  provided  by  the  off-diagonal  matrices,  S  and  L. 
Equation  (1),  which  expresses  a  balance  of  forces,  shows  that  L  relates  the 
pressure  at  the  fluid-structure  interface  to  force  resulting  from  that 
pressure.  In  a  lumped  formulation,  the  terms  of  L  may  be  taken  to  be  the 
interface  surface  areas  associated  with  each  fluid-structure  interface  grid 
point.  In  general,  there  will  be  three  areas  for  each  grid  point  -  one 
for  each  coordinate  direction.  Zienkiewicz^  has  shown  that 

S  =  (6) 

Equation  (5)  then  becomes 


T 

The  L  and  L  matrices  may  be  determined  separately  and  input  directly, 
using  the  UMIG  option  of  NASTRAN.  Thus,  formulation  of  Equation  (7)  (or 
Equation  (5))  is  completed,  and  the  solution  is  carried  out  by  the  program. 
However,  the  presence  of  the  off-diagonal  matrices  destroys  the  symmetry 
of  the  mass  and  stiffness  matrices,  requiring  tlie  use  of  time-consuming 
complex  eigenvalue  routines  in  tlie  solution  of  these  equations. 

The  natural  vibration  frequencies  for  some  flat  plates  with  varying 
levels  of  immersion  were  calculated  by  the  metliod  described  and  comp.ireil 
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with  experimental  results  by  Marcus.^  Good  agreement  was  obtained  when  a 
sufficient  amount  of  water  was  modeled,  and  in  fact  the  results  provide 
a  useful  reference  for  judging  the  size  of  a  fluid  mass  required  to  repre¬ 
sent  an  infinite  body  of  water. 

APPLICATION  TO  PROPELLER  BLADES 

It  was  planned  originally  to  apply  the  analysis  as  represented  by 
Equation  (7)  to  a  supercavitating  propeller  blade,  P3604,  for  which  data 
for  a  finite  element  idealization  were  available  from  a  previous  stress 
study.  The  idealization  consisted  of  155  grid  points  and  16  quadratic 
isoparametric,  twenty-noded ,  three-dimensional  elements  (NASTRAN  data  card 
CIHEX2).  The  purpose  of  this  sample  calculation  was  to  verify  that  appli¬ 
cation  of  this  procedure  to  the  more  complicated  propeller  blade  would  be 
feasible.  It  was  apparent,  however,  that  the  cost  of  such  a  computer  run 
could  be  high.  A  fluid  idealization  was  constructed  of  only  three  layers 
of  elements  and  extended  one  element  beyond  the  edges  of  the  blade.  The 
fluid  idealization,  which  was  thought  to  be  just  large  enough  to  represent 
an  infinite  quantity,  but  only  on  one  side  of  the  blade,  contained  90 
three-dimensional  elements  and  578  grid  points.  A  total  of  263  uncon¬ 
strained  degrees  of  freedom  resulted,  and  tliis  number  would  be  about 
doubled  for  fluid  on  both  sides  of  the  blade.  The  blade  itself  was  con¬ 
sidered  fixed  along  its  base,  so  that  a  total  of  396  degrees  of  freedom 
(3  for  each  unconstrained  grid  point)  remained.  The  total  number  of 
degrees  of  freedom  for  fluid  and  structure  thus  led  to  mass  and  stiffness 
matrices  for  Equation  (7)  of  order  659.  Eigenvalue  extraction  of  these 
large  unsymmetrical  matrices  would  have  required  a  time-consuming  complex 


extraction  routine.  It  was  decided,  therefore,  to  attempt  to  reduce  the 


computer  effort  required  without  any  significant  sacrifice  of  accuracy. 
In  Equations  (2)  and  (7),  the  Q  matrices,  which  contain  a  factor  of 


1/c,  represent  the  effect  of  compressibility  of  the  fluid.  This  effect  is 
not  considered  significant  except  for  high  frequencies  or  in  confined 
fluid  spaces.  Most  analyses,  including  that  of  Everstine,^  have  incor¬ 
porated  the  simplifying  assumption  of  incompressibility;  as  a  result,  1/c, 
and  hence  Q,  are  zero.  This  assumption  was  not  included  in  the  formulation 
of  Equation  (7),  which  was  set  up  for  ready  computer  application,  because 
there  was  little  advantage  in  doing  so. 

IVhen  Q  is  taken  to  be  zero.  Equation  (7)  is  readily  reduced  to 

[M+pLH~^lJ] {u}  +  [K](u}  =  {F;  (g) 

One  additional  significant  operation,  that  of  obtaining  the  inverse 

of  H,  is  required  to  set  up  Equation  (8).  However,  this  and  the  eigenvalue 

extraction  operations  are  performed  on  matrices  about  one-half  the  size 

of  those  of  Equation  (7).  Furthermore,  the  matrices  of  Equation  (8)  are 
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symmetrical.  The  term  pLH  L  is  interpreted  as  the  added  effective  mass 
of  the  fluid.  It  is  clear  that  fluid  on  only  one  side  of  the  propeller 
blade  may  be  idealized,  and  the  added  mass  term  can  be  made  to  represent 
fluid  on  both  surfaces  by  using  a  value  for  fluid  density,  p,  that  is 
double  the  actual  value. 

No  standard  NASTRAN  solution  option,  or  "rigid  format,"  is  available 
to  formulate  and  solve  Equation  (8).  Executive  modules  can  be  manipulated, 
however,  though  the  Direct  Matrix  Abstraction  Program  (DMAP)  capability, 
and  existing  rigid  formats  may  be  altered  in  the  same  way.  These  capa¬ 
bilities  were  used  to  devise  a  procedure  for  implementing  Equation  (8). 


First,  the  fluid  matrix  (or  rather,  (l/p)H)  was  computed  and  made  available 
for  later  reference.  This  fluid  "stiffness"  matrix  was  produced  by  making 
a  static  stress  run  on  the  fluid  mass  with  material  properties  and  boundary 
conditions  as  discussed  above.  ASETl  cards  were  used  to  specify  that  the 
solution  set  be  retained  only  at  the  interface  grid  points,  and  an  "OUTPUTl" 
DMAP  ALTER  instruction  to  specify  that  the  resulting  condensed  stiffness 
matrix  (or  H)  be  saved  on  a  permanent  file.  Another  ALTER  instruction  was 
used  to  stop  further  execution  of  the  problem.  This  run  can  also  be  used 
to  obtain  terms  of  the  L  matrix.  It  was  suggested  previously  that  areas 
associated  with  each  structural  interface  grid  point  can  be  used  as  L 
matrix  terms .  Calculating  these  areas  and  obtaining  their  projections  in 
the  three  coordinate  directions  for  a  propeller  surface  can  be  quite 
difficult.  It  is  much  easier,  and  equivalent,  to  specify  a  unit  pressure 
loading  on  the  interface  and  request  a  printout  of  the  three  equivalent 
nodal  forces  on  these  grid  points.  The  forces  will  thus  be  calculated  in 
a  consistent  manner.  These  forces  can  be  used  to  make  up  the  L  matrix. 

They  are  functions  of  only  the  pressurized  surface  shape  and  the  location 
of  the  grid  points.  Therefore,  since  these  two  quantities  are  identical 
for  both  the  structure  and  the  fluid  at  their  interface,  the  unit  pressure 
load  can  be  applied  to  the  idealized  fluid  interface  in  the  same  run  from 
which  the  H  matrix  is  obtained.  The  equivalent  nodal  forces  can  then  be 
entered  on  DMIG  cards  as  L  matrix  terms.  Since  the  fluid  run  has  only  one 
degree  of  freedom,  while  all  three  translational  degrees  are  needed  to 
obtain  the  three  components  for  each  surface  grid  point,  the  equivalent 
force  matrix  must  be  printed  out  before  constraints  are  applied.  This 
printout  is  requested  through  another  ALTER  instruction. 
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To  summarize  the  procedure:  a  static  stress  run  is  started  on  the 
fluid  mass  with  unit  pressure  applied  to  the  interface.  Ail  interface  grid 
points  are  specified  on  ASETl  cards.  The  stiffness  matrix  (after  con¬ 
straints  and  condensation  are  applied)  is  placed  on  a  permanent  file  via 
the  required  control  card.  An  OUTPUTl  ALTER  instruction  and  the  equivalent 
nodal  forces  are  then  printed  out,  at  which  point  the  calculation  is 
stopped.  The  matrix  condensation  is  not  theoretically  required,  but  is 
done  to  save  space;  it  does  not  affect  the  final  results. 

Equation  (8)  may  be  set  up  for  the  propeller  and  solved  to  obtain  ei¬ 
ther  the  frequency  response  or  the  dynamic  response  to  applied  time-depen¬ 
dent  forces  by  running  the  propeller  data  with  the  additional  data  already  • 
obtained.  The  natural  frequencies  are  obtained  by  specifying  a  modified 
natural  modes  calculation  with  F  null  (Rigid  Format  3);  the  dynamic 
response  is  obtained  by  specifying  a  transient  response  calculation  with 
F  input  (Rigid  Format  9).  Each  of  these  rigid  formats  requires  its  own 
ALTER  routine,  although  each  routine  works  in  essentially  the  same  way, 
and  the  additional  input  is  handled  similarly.  Structural  grid  point  data 
must  include  one  extra  scalar  point  (SPOINT)  for  each  fluid  intertaLC  grid 
point  listed  on  ASETl  cards  and  thus  retained  as  a  degree  of  freedom  in 
the  H  matrix  calculation.  These  scalar  points  are  dummy  grid  points  which 
have  no  coordinates,  but  have  one  allowable  degree  of  freedom  each  (the 
zero  component).  They  are  established  by  listing  their  ident i t icat ion 
numbers  on  SP01^T  cards.  These  identification  numbers  are  (for  conven¬ 
ience)  the  fluid  interface  grid  point  numbers,  and  they  must  be  dillerent 
from  ail  propeller  grid  point  numbers.  Tlie  purpose  of  the  scalar  points 
is  to  expand  the  dimensions  of  the  structural  stiffness  matrix,  belore 
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T 

constraints  are  applied,  sufficiently  to  accept  the  L  (or  L  )  matrix  which 
is  also  input  to  the  stiffness  matrix  via  DMIG  cards.  It  was  found  more 
advantageous  to  require  the  L  matrix;  the  L  matrix  can  readily  be  derived 
from  it  by  transposition.  The  terms  of  the  matrix  are  specified  on  DMIG 
cards  ,  each  of  which  relates  one  of  the  three  nodal  force  components  per 
unit  pressure  at  a  wetted  propeller  grid  point,  obtained  from  the  fluid 

calculation  with  the  zero  component  of  the  corresonding  scalar  point.  If 

T 

the  algebraic  signs  of  all  of  the  L  terms  are  reversed,  those  of  L  will 

—  1  T 

also  be  reversed,  and  the  added  mass  term,  pLH  L  ,  will  be  unchanged; 

T 

therefore  a  sign  error  for  the  L  matrix  as  a  whole  has  no  effect.  Since 

only  the  interface  fluid  points  are  input  here  (as  scalar  points),  the  zero 

rows  of  L  are  eliminated  so  that  the  matrix  conforms  with  the  condensed 
_2. 

version  of  H  .  The  (l/p)H  matrix  is  obtained  from  the  file  on  which  it  has 
been  stored  by  including  the  required  control  card  and  an  INPUTTl  ALTER 
instruction. 

Two  partitioning  matrices,  named  PVECl  and  PVECDl ,  must  also  be  input 

via  DMI  cards.  These  are  column  matrices  having  the  dimensions  (N+S)xl.  K 

is  the  total  number  of  structural  degrees  of  freedom  (number  of  grid  points 

x6)  and  S  is  the  number  of  scalar  points.  For  PVECl,  the  first  N  terms 

are  zero,  and  the  remaining  S  terms  are  l.U.  PVECl  is  used  to  partition 

T 

the  stiffness  matrix  to  obtain  L  ,  used  with  (l/p)H  by  the  ALTER  routine  to 

-1  T 

form  the  added  mass  term  pLH  L  ,  which  is  then  added  to  the  structural 
mass  matrix.  Note  that  the  inverse  of  (l/p)H  is  pH  For  the  PVEDCl  matrix, 
all  (N+S)  terms  are  zero.  PVECDl  is  used  to  form  a  "dummy"  partition  on 
tlie  "augmented"  mass  matrix,  which  now  includes  the  added  mass  terms,  so 
that  the  program  will  correctly  perceive  this  matrix  as  symmetrical. 
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;)ne  tinal  type  of  additional  input  is  required:  all  scalar  points  must 
be  listed  on  single  point  constraint  (SPC)  cards.  These  constraints  are  in 
addition  to  whatever  constraints  are  applied  to  the  propeller.  It  is  sug¬ 
gested  that  2U-node  three-dimensional  C1HEX2  elements,  with  two  integration 
points  specified,  be  used  to  represent  both  the  fluid  and  the  propeller, 
and  that  tlie  consistent  mass  matrix  formulation  be  specified  for  the  pro¬ 
peller. 

The  procedure  just  described  can  be  used  to  obtain  vibration  modes 
and  dynamic  response  of  propellers  under  a  variety  of  conditions  of  sub¬ 
mergence.  Partial  submergence  at  any  given  instant  can  be  represented 
during  a  fluid  run  by  listing  on  ASETl  cards  only  the  fluid  grid  points 
wliich  correspond  to  submerged  points  on  the  propeller.  During  a  structural 
run  these  points  are  listed  on  SPOINT  cards.  Also,  the  dimensions  of  the 
PVECi  and  PVECDl  column  matrices  must  be  adjusted  by  making  the  number  of 
terms  with  value  of  1.0  at  the  end  of  the  PVECI  column  equal  to  the 
number  of  wetted  grid  points.  Wetting  of  both  surfaces  of  the  blade  is 
simulated  by  using  double  the  actual  value  of  fluid  dcusit;,,  .  ,  for  the 
fluid  run.  Wetting  of  only  one  surface,  sucli  as  occurs  during  supercavi¬ 
tation,  is  represented  by  use  of  the  actual  value  of  c. 

This  procedure  can  be  extended,  witliout  further  DM/VP  program  modifi¬ 
cations,  to  analyze  dynamic  response  of  partially  submerged  propellers 
and  to  take  into  account  the  progressively  changing  degree  of  submergence. 
Tlie  extended  procedure  consists  essentially  of  running  the  analysis  in 
separate  steps,  one  for  each  arbitrarily  cliosen  stage  of  submergence.  A 
convenient  way  to  determine  a  new  stage  ol  submergence  is  to  consider  the 
wetting  ot  a  new  row  of  propeller  elements.  Tliis  extended  procedure  was 
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devised  and  was  applied  to  an  existing  propeller  for  which  some  experi¬ 
mental  and  analytical  data  were  available.  The  details  of  this  application 
will  be  presented  later.  The  finite  element  idealization  of  this  blade 
was  laid  out  so  that  one  series  of  element  boundaries  coincided  with  water 
lines  on  the  propeller  blades  as  it  progressed  into  the  water  during  each 
revolution.  There  were  two  such  boundaries  in  this  rather  coarse  ideali¬ 
zation,  defining  three  rows  of  elements;  hence,  three  stages  of  submergence 
were  assumed. 

The  fluid  run,  to  obtain  the  condensed  (l/p)H  matrix,  was  run  once  for 
each  stage  of  submergence.  For  each  run,  the  solution  set  of  wetted  grid 
points,  as  defined  by  AbETl  cards,  contained  only  the  points  presently 
wetted,  i.e.,  points  just  submerged  plus  all  points  previously  submerged. 
Similarly,  unit  pressure  was  applied  to  the  corresponding  wetted  element 
faces  to  obtain  the  L  matrix  terms  for  the  same  wetted  points.  Finally, 
those  interface  grid  points  not  yet  submerged  had  to  be  considered  part 
of  the  fluid  boundary  and  were  constrained.  Constraints  were,  of  course, 
removed  from  those  points  as  they  became  submerged  in  subsequent  stages  and 
added  to  those  included  on  ASETl  cards.  Each  of  the  (l/p)H  matrices  thus 
produced  and  stored  on  a  permanent  file  was  assigned  a  separate  name;  in 
this  case  the  names  used  were  HMATRIXSTEP 1 ,  HMATR1XSTEP2 ,  and  HMATRIXSTEP3 . 
Sucli  separate  identification  of  the  L  matrices  was  not  necessary,  as 
they  are  manually  punched  on  DMIG  cards  and  input  separately  to  the  struc¬ 
tural  computer  runs.  The  structural  run  in  this  case  was  also  carried 
out  once  for  each  stage  of  submergence,  but  all  except  the  first  were 
restart  runs.  Thus,  ail  but  the  last  run  was  checkpointed. 
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The  series  of  runs  is  set  up  as  follows.  First,  time-dependent 
forces  are  specified  for  the  selected  grid  points  by  means  of  the 
required  DAREA,  DLOAD,  TABLEDl ,  and  TLOADi  cards.  The  time  spanned  by 
these  forces  covers  the  entire  time  for  which  the  analysis  is  applied. 

The  time  at  which  the  force  is  initiated  at  each  selected  grid  point  is 
dictated  by  the  speed  of  rotation  of  the  propeller  and  by  the  location  of 
the  point.  The  user  must  also  specify  the  length  of  each  integration  time 
step,  and  the  number  of  steps,  on  a  TSTEP  card.  The  calculation  is  ter¬ 
minated  when  the  time  defined  by  the  product  of  time  step  x  number  of  steps 
on  the  TSTEP  card  has  elapsed.  If  a  uniform  time  step  is  used  throughout, 
the  number  of  steps  may  be  specified  to  make  each  run  cover  just  the  length 
of  time  required  for  the  water  line  to  progress  from  the  first  to  the 
second  boundary  defining  the  present  stage  of  submergence.  The  first  run 
and  each  restart  then  require  a  different  TSTEP  card  (unless,  by  coinci¬ 
dence,  each  stage  of  submergence  lasts  the  same  length  of  time).  Other 
changes  must  also  be  made  between  runs.  Scalar  points  for  all  interface 
fluid  grid  points  which  are  to  be  wetted  at  any  time  during  the  analysis 
must  be  input  (on  SPOINT  cards)  for  the  first  run,  and  this  does  not  change 
for  succeeding  restarts.  The  same  is  true  of  the  scalar  point  constraints. 
The  PVECOl  matrix  is  also  input  as  before,  with  column  dimension  equal  to 
six  times  the  number  of  propeller  grid  points  plus  the  number  of  scalar 
points,  and  it  is  not  changed.  The  PVECl  column  matrix  is  now  replaced 
by  a  series  of  matrices,  one  for  each  structural  run,  and  they  may  be 
designated  PVECl,  PVEC2 ,  and  PVEC3 .  The  column  dimension  is  the  same  for 
all  of  these,  and  is  the  same  as  that  of  PVECDl .  The  number  and  location 
ot  the  "l.U"  terms,  however,  must  correspond  to  the  toral  set  of  wetted 
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devised  and  was  applied  to  an  existing  propeller  for  which  some  experi¬ 
mental  and  analytical  data  were  available.  The  details  of  this  application 
will  be  presented  later.  The  finite  element  idealization  of  this  blade 
was  laid  out  so  that  one  series  of  element  boundaries  coincided  wit)i  water 
lines  on  the  propeller  blades  as  it  progressed  into  the  water  during  each 
revolution.  There  were  two  such  boundaries  in  tliis  rather  coarse  ideali¬ 
zation,  defining  three  rows  of  elements;  hence,  three  stages  of  submergence 
were  assumed. 

The  fluid  run,  to  obtain  the  condensed  (l/p)H  matrix,  was  run  once  for 
each  stage  of  submergence.  For  each  run,  the  solution  set  of  wetted  grid 
points,  as  defined  by  ASETl  cards,  contained  only  the  points  presently 
wetted,  i.e.,  points  just  submerged  plus  all  points  previously  submerged. 
Similarly,  unit  pressure  was  applied  to  the  corresponding  wetted  element 
faces  to  obtain  the  L  matrix  terms  for  the  same  wetted  points.  Finally, 
those  interface  grid  points  not  yet  submerged  had  to  be  considered  part 
of  the  fluid  boundary  and  were  constrained.  Constraints  were,  of  course, 
removed  from  those  points  as  they  became  submerged  in  subsequent  stages  and 
added  to  those  included  on  ASETl  cards.  Each  of  the  (1/())H  matrices  thus 
produced  and  stored  on  a  permanent  tile  was  assigned  a  separate  name;  in 
this  case  the  names  used  were  UMATRIXSTEPI ,  HMATK1XSTEF2 ,  and  HMATRIXSTEF • . 
Such  separate  identification  of  tlie  L  matrices  was  not  necessary,  as 
they  are  manually  punched  on  DMIC  cards  and  input  separately  to  the  struc¬ 
tural  computer  runs.  The  structural  run  in  this  case  was  also  carried 
out  once  for  each  stage  of  submergence,  but  all  except  the  first  were 
restart  runs.  Thus,  all  but  the  last  run  was  check[)ointed . 


The  series  of  runs  is  set  up  as  follows.  First,  time-dependent 
forces  are  specified  for  the  selected  grid  points  by  means  of  the 
required  DAREA,  DLOAD,  TABLEDl ,  and  TLOADi  cards.  The  time  spanned  by 
these  forces  covers  the  entire  time  for  which  the  analysis  is  applied. 

The  time  at  which  the  force  is  initiated  at  each  selected  grid  point  is 
dictated  by  the  speed  of  rotation  of  the  propeller  and  by  the  location  of 
the  point.  The  user  must  also  specify  the  length  of  each  integration  time 
step,  and  the  number  of  steps,  on  a  TSTEP  card.  The  calculation  is  ter¬ 
minated  when  the  time  defined  by  the  product  of  time  step  x  number  of  steps 
on  the  TSTEP  card  has  elapsed.  If  a  uniform  time  step  is  used  throughout, 
the  nuiaber  of  steps  may  be  specified  to  make  each  run  cover  just  the  length 
of  time  required  for  the  water  line  to  progress  from  the  first  to  the 
second  boundary  defining  the  present  stage  of  submergence.  The  first  run 
and  each  restart  then  require  a  different  TSTEP  card  (unless,  by  coinci¬ 
dence,  each  stage  of  submergence  lasts  the  same  length  of  time).  OtJier 
changes  must  also  be  made  between  runs.  Scalar  points  for  all  interface 
tiuid  grid  points  which  are  to  be  wetted  at  any  time  during  the  analysis 
must  be  input  (on  SPUlNT  cards)  for  the  first  run,  and  this  does  not  change 
for  succeeding  restarts.  The  same  is  true  of  the  scalar  point  constraints. 
The  PVECDl  matrix  is  also  input  as  before,  with  column  dimension  equal  to 
six  times  the  number  of  propeller  grid  points  plus  the  number  of  scalar 
points,  and  it  is  not  changed.  The  PVECl  column  matrix  is  now  replaced 
by  a  series  of  matrices  ,  one  for  each  structural  run,  and  they  may  be 
designated  PVECl,  PVEC2 ,  and  PVEC3.  The  column  dimension  is  the  same  for 
dll  of  these,  and  is  the  same  as  that  of  PVECDl.  The  number  and  location 
of  the  "l.u"  terms,  however,  must  correspond  to  the  total  set  of  wetted 


interface  points  as  represented  by  the  scalar  points  during  each  run.  The 
DMAP  ALTER  instruction  which  utilized  FVEC  must  also  be  changed  from  run 
to  run  to  call  for  the  currently  available  PVECi  matrix.  Another  ALTER 
instruction,  as  well  as  a  control  card,  must  be  changed  tor  each  run  to 
obtain  the  required  (l/p)H  matrix  for  each  stage.  Finally,  each  run  must 
also  have  a  new  L^  matrix  input  via  DMKJ  cards.  Since  tliis  information 
is  part  of  the  bulk  data,  the  previous  L^  matrix  must  be  eliminated  during 
eacli  restart. 

Sample  input  for  a  vibration  problem  and  for  a  transient  response 
problem  are  presented  in  tlie  appendix. 

VIBRATION  ANALYSIS  OF  P3604  PROPELLER 
The  vibration  analysis  procedure  was  first  applied  to  a  supercavitating 
propeller,  P3b04,  for  which  finite  element  idealization  data  were  already 
available  from  a  previous  stress  study.  The  blade  was  divided  into  lb 
isoparametric,  20-noded ,  three-dimensional  elements,  with  135  grid  points 
(Figure  1).  A  quantity  of  water  tliree  elements  thick,  extending  one  ele¬ 
ment  beyond  the  edges  of  the  blade,  and  generally  following  t  lie  twist  ol 
the  blade,  was  assumed.  The  elements  were  roughly  cubical,  and  the  thick¬ 
ness  of  the  water  modeled  was  about  one-half  to  two-thirds  the  length  ot 
the  blade.  This  size  was  judged  to  be  sufficient  to  represent  an  infinite 
body  of  water,  based  on  a  study  of  the  results  presented  by  Marcus.'^  With 
the  procedure  used  here,  larger  representations  could  be  investigated  at 
fairly  moderate  cost.  The  double  value  of  water  density  was  input  to 
represent  water  on  both  sides  of  the  blade.  This  fluid  idealization  con¬ 
sisted  of  578  grid  points  and  90  elements.  A  total  ot  2()3  degrees  of 
freedom  remained  after  constraints  were  applied.  The  fluid-structure 
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ioterface  surface  contained  65  grid  points,  which  were  retained  in  the 


T 

H  matrix.  The  L  matrix  terms  numbered  65x3  or  195.  However,  nine  of 

the  interface  points  were  adjacent  to  the  base  grid  points  of  the  propeller, 

T 

wliich  were  later  to  be  constrained.  The  corresponding  9x3,  or  27  L  matrix 

terms  were  therefore  not  Included  on  DMIG  cards.  Space  was  later  reserved 

for  them  in  the  propeller  matrices,  however,  as  65  scalar  points  were 
T 

input,  but  the  L  boundary  terms  wliich  were  not  input  were  implicitly  zero. 
It  was  not  necessary  to  leave  out  these  terms;  they  would  be  eliminated 
later  by  constraints.  With  this  procedure,  the  pressures  in  the  water  at 
the  nine  base  grid  points  were  not  constrained  to  zero,  even  though  deflec¬ 
tions  of  the  corresponding  propeller  grid  points  were  constrained.  Tlie 
tive  lowest  natural  frequencies  of  the  blade  in  water  were  calculated  and 
are  presented,  with  the  corresonding  calculated  in-air  frequencies,  in 
Table  1.  The  modal  deflection  data  indicated  that  the  first  mode  of 
vibration  consisted  of  nearly  pure  spanwise  bending,  while  the  other  modes 
contained  a  combination  of  torsion  and  chordwise  bending.  Although  no 
experimental  data  exist  to  compare  with  this  sample  calculation,  the 
results  appear  to  be  consistent  with  experimental  results  for  more  conven¬ 
tional  blades,  particularly  for  the  first  mode.  It  seems  rather  remarkable 
that  the  ratios  of  frequencies  in  water  to  those  in  air  are  fairly  constant, 
ranging  from  U.46  to  0.55,  for  all  modes.  This  was  not  expected  because 
of  the  large  differences  in  thickness  of  the  blade  and  the  variety  ul 
possible  vibration  modes.  This  phenomenon  may  not  be  characteristic  of 
other  blade  geometries. 

The  version  of  NASTRAN  used  for  tliese  analyses  was  Level  17.  I’roce- 
Uures  for  both  the  vibration  and  transient  analyses  are  listed  in  the 
Appendix. 
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TRANSIENT  RESPONSE  ANALYSIS  OE  SES-IUUB  PROPELLER 


After  this  work  was  underway,  some  fuii  scale  test  results,  includinj; 
both  vibration  and  transient  response  data,  became  available^’^  tor  tlie  SES- 
lOOB  controllable  pitch,  partially  submerged,  supercavitatlng  propeller 
(number  7016.3).  Allison's  work^  also  contains  finite  element  analyses  of 
the  vibration  modes  and  transient  response  of  this  blade,  utilizing  27b 
grid  points  and  379  plate  bending  elements.  Calculated  transient  forces 
were  used,  and  a  constant  pliysical  mass  ot  water  ol  approxiiiMtely  conical 
shape  circumscribing  the  blade,  was  taken  as  added  effective  mass.  One-lialf 
of  the  imaginary  cone  of  water  was  actually  used  to  simulate  full  cavita¬ 
tion.  Transient  response,  in  the  form  of  measured  stresses  at  a  gage  near 
the  root  of  the  blade,  was  plotted  versus  time  lor  the  i)eriod  during  which 
the  blade  enters  the  water  and  reaches  maximum  submergence. 

The  SES-lUt)B  propeller  was  selected  for  this  transient  response  analy¬ 
sis  because  of  the  availability  of  geometrical  .ind  test  data.  Due  to  the 
developmental  nature  of  this  ap|)licat  itui,  however,  and  to  time  and  cost 
constraints,  the  idealization  constructed  was  rather  coarse,  consisting  of 
96  grid  points  and  nine  C1HEX3  isopar.imet r ic  three-dimensional  elements. 

A  detailed  study  of  this  blade  would  reijuire  a  finer  idealization,  particu¬ 
larly  so  because  it  contained  an  ,ninex,  or  d  iscoitt  inui  t  y ,  in  the  pressure 
surface  curvature.  The  exercise  ol  engineering  judgment  , it  tempted  to  com¬ 
pensate  for  tills  coarseness  in  constructing  the  idealization,  in  which, 
for  instance,  thicknesses  at  grid  points  do  not  all  coincide  exactly  with 
those  on  the  actual  blade.  tine  ol  the  constraints  on  placement  of  element 
boundaries  was  to  require  that  one  set  ol  boundaries  coincide  with  water 


lines  at  successive  assumed  stages  of  submergence  of  the  blade.  This 


constraint  would  not  be  necessary  with  a  fine  idealization. 

A  view  of  the  blade,  looking  in  the  forward  direction,  is  siiown  in 
Figure  1.  A  series  of  water  lines  at  various  stages  of  submergence  was 
constructed  from  data  reported  by  Lewis^  and  drawn  on  Figure  1.  An  ideali¬ 
zation  was  constructed  with  two  of  the  element  boundaries  approximately 
coinciding  with  two  water  lines,  and  this  is  shown  in  Figures  2  and  3.  In 
these  figures  the  water  lines  are  shown  slightly  curved  because  the  view 
of  the  propeller  is  flattened.  The  location  of  the  strain  gage  from  which 
Allison's  transient  response  data^  were  obtained  is  also  shown  in  Figure  2. 
For  the  first  arbitrarily-defined  stage  of  submergence,  elements  7  through 
9  are  submerged;  for  the  second  stage,  elements  4  through  9  are  submerged; 
and  for  the  third  stage,  all  elements  are  submerged.  The  fluid  mass 
idealized  for  this  case  was  of  approximately  the  same  size  and  shape  rela¬ 
tive  to  the  blade  as  for  the  P3604  vibration  analysis,  but  the  fluid  for 
the  SFS-IOOB  was  divided  into  only  two  layers.  An  outline  of  one  face  of 
the  fluid  mass,  showing  the  fluid-structure  interface,  is  sliown  in  Figure  4. 
The  entire  fluid  nujss  is  shown  in  Figure  5. 

The  operating  conditions  for  which  tills  calculation  was  curried  out 
were  the  design  conditions  of  83  knots,  1870  rpm,  and  33-degree  blade 
pitch.  Blade  pressures  are  discussed  by  Allison’  and  l.ewis,*^’  both  ot  whom 
considered  the  pressure  at  each  point  to  be  applied  in  three  stages.  first 
was  an  impact  pressure,  compvited  according  to  a  theory  by  Pierson,^  fol¬ 
lowed  by  an  entry  pressure  predicted  from  a  tiieory  by  Wang.^  Final  iy,  a 
steady  value  was  reached,  tlie  value  predicted  from  a  tiieory  by  Yim.  The 
same  procedure  was  followeii  lien-  to  calculate  forces  to  applv  to  grid 


points  73,  75,  A7 ,  A9,  21,  and  23.  These  are  mid-side  nodes  and  are  so 
situated  that  force  is  applied  to  each  when  tlie  propeller  reaches  the  mid¬ 
point  of  its  stage  of  submergence.  The  pressure  versus  time  plot  at  each 
of  these  grid  points  was  derived  from  data  given  by  Lewis. ^  These  plots 
are  presented  in  Figure  6  and  show  that  these  are  essentially  pulse  pres¬ 
sures  with  no  oscillating  component. 

The  transient  response  analysis  consisted  ol  one  fluid  and  one  pro¬ 
peller  computer  run  for  each  stage  of  submergence,  a  total  of  six  runs. 
Stresses  are  printed  out  in  NASTKAN  at  the  grid  points  and  at  the  element 
centroids.  Unfortunately,  it  was  not  feasible  to  place  a  grid  point 
exactly  where  experimental  data  were  available.  Therefore,  stresses 
computed  at  three  grid  points  (22,  36,  and  39)  which  surround  this  point 
(stiown  as  Gage  1  in  Figure  2)  were  plotted  and  the  plots  are  shown  in 
r'igurcs  7,  g,  and  9.  Static  stresses  were  also  calculated  and  plotted  on 
tliese  same  figures.  Tlie  three  portions  of  the  dynamic  forces  were  averaged 
and  applied  at  tlie  same  points  to  obtain  these  static  stresses.  Tlie  dyna¬ 
mic  pressures  were  also  applied,  but  with  the  steady  state  value  truncated, 
to  the  propeller  in  air,  and  the  results  are  given  in  Figures  lU,  11,  and 
12.  The  plot  of  dynamic  stresses,  reproduced  from  Allison's  paper, is 
presented  here  as  Figure  13.  Natural  vibration  frequencies  for  the 
SES-IOUB  propeller  at  the  three  stages  of  submergence  were  also  calculated 
and  are  presented  in  Table  2. 

DISCUSSION  OF  SES-IUOB  ANALYSIS  RESULTS 

The  primary  purpose  of  ti)is  analysis  was  to  verify  the  computational 
procedures  developed  and  to  produce  results  that  appear  reasonable  and  that 
provide  some  insiglit  into  the  beliavior  ol  svicli  propellers.  It  was  also 
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desired  that,  despite  the  coarseness  of  the  idealization,  some  meaningful 


comparisons  could  be  made  witii  experimental  data. 

The  plots  of  transient  stresses  in  the  propeller  in  water  (Figures  7, 
8,  and  9)  show  that  the  peak  stresses  at  the  three  locations  plotted  range 
from  21500  to  32000  psl.  The  measured  peak  stress  shown  in  Figure  13  is 
about  31000  psi.  This  apparent  close  agreement  is  probably  fortuitous. 

The  stress  amplitudes  shown  in  the  analytical  plots  vary  with  a  predominant 
frequency  equal  to  the  lowest  natural  frequency  of  the  propeller,  as  veri¬ 
fied  by  the  natural  frequency  analysis.  What  appears  to  be  a  similar  pre¬ 
dominant  vibrational  component,  with  nearly  the  same  frequency,  is  dis¬ 
played  in  the  experimental  plot  (Figure  13).  Substantial  amplitude  damping 
is  evident  in  the  experimental  plot,  however,  even  within  each  complete 
period  of  vibration.  Both  the  analytical  and  experimental  plots  indicate 
that  the  blade  vibrates  through  several  complete  cycles  before  becoming 
fully  submerged.  Thus,  pressure  due  to  submergence  in  the  later  stages 
are  applied  when  the  existing  calculated  amplitudes  of  deflection  are  in 
error,  due  to  the  lack  of  damping.  Furthermore,  the  calculated  phase  could 
be  in  error  due  to  the  probable  slight  inaccuracy  of  the  calculated  natural 
frequencies.  Calculated  response  due  to  these  later  water  entry  pressures 
may  therefore  detract  from,  rather  than  improve,  accuracy  of  the  overall 
dynamic  response.  The  calculated  stress-time  plots  show  no  apparent  reac¬ 
tion  to  these  later  stage  pressures,  however.  This  is  evidently  because 
the  last  part  of  the  propeller  to  enter  tlie  water  is  the  portion  close  to 
the  liLib  where  lower  pressures  are  generated  at  iiiqiact  (due  to  the  smaller 

*Damping  was  not  accounted  for  in  ttio  analysis,  hut  could  be,  as  discussed 
later  . 
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radii  and  hence  lower  linear  velocity  resulting  from  rotation)  and  impact 
pressures  produce  smaller  bending  moments  due  to  the  smaller  effective 
moment  arms.  The  plot  of  calculated  stress  versus  time  is  probably 
meaningful  quantitatively  only  to  the  first  peak  if  damping  is  not  ac¬ 
counted  for.  Fortunately,  this  portion  of  the  curve,  or  more  particularly 
the  value  of  this  peak,  is  of  primary  interest. 

Two  other  features  of  these  curves  are  significant.  First,  peak 
stress  values  are  practically  the  same  for  the  in-water  and  the  In-air 
cases.  Second,  these  peak  values  are  on  the  order  of  two  to  three  times 
the  static  values.  This  suggests  that,  for  this  approximately  pulse-type 
pressure  application,  the  blade  behaves  somewhat  like  a  weight  suddenly 
suspended  by  a  spring:  the  weight  oscillates  about  the  static  equilibrium 
point  with  an  amplitude  double  the  static  deflection.  However,  the  ratios 
of  the  dynamic  to  static  stresses  in  the  propeller  blade  probably  vary  at 
different  locations  because  of  the  complicated  shape.  For  this  particular 
example,  a  Fourier  series  representation  of  the  assumed  pressure  pulses 
apparently  would  not  contain  any  strong  components  near  any  of  the  natural 
frequencies  in  water  or  in  air.  T'le  damping  effect  of  water  no  doubt 
serves  to  decrease  peak  dynamic  stresses  slightly,  but  the  main  effect  of 
submersion  in  water,  apart  from  transmitting  force  to  the  blade,  appears 
to  be  in  altering  its  natural  frequencies.  If  a  detailed  study  of  a  pro¬ 
peller  were  to  be  made,  it  would  probably  be  wise  to  obtain  stresses  at 
all  grid  points  (at  least  all  the  points  on  one  surface)  because  a  rea- 
sonance  magnification  might  occur  locally,  and  hardly  show  up  at  the  root 
of  the  blade. 
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The  t undamentai ,  or  lowest,  natural  frequencies  of  the  blade  in  air 


and  in  water  may  be  derived  from  the  stress  versus  time  plots  since  the 
time  record  is  probably  long  enough.  All  three  plots  for  the  in-air  case 
display  frequencies  of  about  65U  Hz,  slightly  higher  than  the  calculated 
natural  frequency  of  624  Hz  shown  in  Table  2.  The  plots  indicate  tliat  in 
water  tlie  blade  rings  at  a  fairly  constant  frequency  of  about  510  Hz,  which 
is  close  to  the  natural  frequency  of  5''4  Hz  calculated  for  the  blade  under 
its  first  stage  of  submergence.  During  most  of  the  time  covered  by  the 
plots,  liowever,  tlie  blade  is  fully  submerged,  and  one  would  expect  the  fre¬ 
quency  to  fall  to  about  420  Hz.  A  smaller  initial  stage  of  submersion,  for 
which  tlie  fundamental  frequency  is  higher,  might  produce  a  corresondingly 
higlier  ringing  frequency.  In  this  case,  however,  the  initial  load  would 
also  be  smaller,  being  applied  to  the  smaller  area  of  the  blade,  and  might 
not  predominate  -iS  it  does  in  the  present  case. 

Damping  (either  hydrodynamic  or  structural)  could  readily  be  accounted 
for  in  this  analysis,  but  the  values  of  damping  terms  would  have  to  be  de¬ 
termined  by  tlie  user,  probably  empirically.  Ail  damping  would  be  accounted 
for  in  the  structural  computations.  Since  damping  terms  can  be  derived 
from  the  reJative  amplitudes  of  successive  cycles  of  a  freely  vibrating 
structure,  plots  such  as  Figure  13,  in  wliich  the  attenuation  of  amplitude 
after  initial  response  is  clearly  illustrated,  could  probably  be  used  tor 
this  [lurpose.  The  principal  effect  of  damping  on  peak  calculated  stresses 
would  probably  occur  when  two  or  more  substantial  dynamic  forces  are 
applied,  separated  by  time  intervals  of  tlie  order  of  one-half  a  1  undamentai 
period . 
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It  would  be  pointles*!  to  carry  the  analysis  further  in  time  to  cover 
the  emergence  of  the  blade  from  the  water,  since  the  dynamic  forces  applied 
at  that  time  would  result  essentially  from  releasing  the  relatively  small 
steady  state  submergence  pressures.  These  forces  would  be  much  less  severe 
than  those  at  entry. 

It  has  been  assumed  in  this  analysis  that  dynamic  behavior  of  a  pro¬ 
peller  can  be  adequately  described  by  analysis  of  a  single  blade,  fixed 
along  its  base.  Errors  due  to  this  assumption  may  be  significant  in  some 
cases,  particularly  for  propellers  witliout  relatively  massive  hubs.  The 
errors  could  be  Introduced  in  two  ways.  First,  flexibility  of  the  hub 
violates  the  assumption  of  complete  fixity,  leading  ^o  lower  natural  fre¬ 
quencies  than  calculated.  This  effect  was  noted  by  Allison^  in  connection 
with  the  SES-IUUB  propeller  and  with  tests  of  a  single  blade  mounted  in 
a  test  jig.  Second,  hub  flexibility  and  lack  of  great  mass  may  permit 
a  significant  degree  of  dynamic  interaction  among  the  blades.  It  is  also 
possible  that  in  some  cases  significant  interaction  may  occur  through 
the  water  between  the  blades. 

The  computational  procedures  followed  here  could  no  doubt  be  simpli¬ 
fied  and  automated  to  some  extent  with  the  aid  of  further  executive  module 
programming. 

CONCLUSIONS 

1.  Rational  procedures  have  been  presented  for  calculating  the  natural 
vibration  modes  and  transient  response  of  partially  or  fully  submerged, 
supercavitating  or  conventional  marine  propellers.  Cost  of  these  analyses 
should  be  moderate,  and  they  require  only  use  of  a  standard  NASTRAN  program. 

2.  Both  the  vibration  modes  and  the  transient  response  analyses,  particu¬ 
larly  when  progressive  submergence  is  treated,  need  more  verification. 


3.  For  the  essentially  pulse-shaped  press ure-time  profile,  as  used  here  in 
the  SES-IOOB  propeller  example,  peak  root  stresses  were  predicted  equally 
well  whether  the  blade  was  considered  in  air  or  in  water.  This  may  not 
have  been  true  of  stresses  elsewhere  on  the  blade.  Also,  it  may  have  been 
by  chance  that  no  resonances  of  the  blade  were  excited  by  components  of 
the  assumed  pressures.  If  time-dependent  pressures  with  oscillatory  com¬ 
ponents  were  applied,  it  is  expected  that  peak  stresses  would  be  affected 
by  submergence  in  water. 

A.  Calculated  stresses  indicated,  as  did  published  experimental  results, 
that  blade  root  stresses  due  to  pulse  pressures  can  be  two  or  more  times  as 
high  as  calculated  stresses  based  on  the  same  pressures  applied  statically, 
b.  Detailed  studies  of  stresses  and  deflections  throughout  typical  pro¬ 
peller  blades,  utilizing  finer  finite  element  meshes  would  be  useful  in 
providing  more  information  about  response  of  typical  propellers  under 
various  conditions. 

6.  When  damping  was  not  taken  into  account,  transient  response  after  about 
ttie  first  half  of  a  fundamental  vibration  period  was  quite  unrealistic, 

but  peak  stress  values  were  apparently  not  affected.  Damping  constants 
need  to  be  established,  perhaps  through  an  empirical  approach.  Tlie 
apparent  failure  of  the  blade  ringing  frequency  to  fall  promptly  to  the 
fully  submerged  value  needs  further  investigation. 

7.  It  would  be  worthwhile  to  simplify  and  automate,  whenever  possible, 
the  procedures  used  here. 
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Figure  6  -  Time-Dependent  Blade  Pressures 
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Figure  10  -  Calculated  Radial  Stress  at  Grid  Point  22  of 
Propeller  Blade  in  Air  with  Applied  Dynamic  Forces 
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TABLE  1  -  CALCULATED  NATURAL  FREQUENCIES  OF  P3604  PROPELLER  BIADE 


Frequency 


Vibration 
Mode,  N 


Calculated 

Frequenc ies 

In  Air 

In  Water 

461 

225 

1228 

671 

1680 

774 

2397 

1171 

2613 

1315 

„  .  in  Water 

Ratio,  — — -r: - 

sn  Air 


TABLE  2  -  NATURAL  FREQUENCIES  OF  SES-IOOB  PROPELLER  BLADE 


V ibrat ion 
Mode,  N 

Calculated 

Frequenc ies 

1  n  Air 

In  Water 

Stage  1 

Stage  2 

Stage  2 

1 

624 

524 

425 

420 

2 

1399 

1037 

924 

916 

3 

1679 

1316 

1206 

1158 

2488 

1893 

1  581 

5 

2612 

204  5 

1815 

1797 

APPENDIX  -  INPUT  DATA  PPEPAKATiUN 


Input  required  for  the  vibration  and  transient  response  analyses 
in  this  report  will  be  described  using  ttie  t>ES-iUUb  calculations  as 
examples.  Inputs  described  will  be  those  required  for  the  fluid  II  and 
L  matrix  runs,  propeller  vibration  runs  in  water,  and  propeller  transient 
response  runs,  with  progressive  submer"ence  of  the  blade.  The  fluid  matrix 
is  the  same  for  either  the  vibration  or  the  transient  response  runs. 

NASTRAN  input  decks  consist  of  control  cards,  an  Executive  Control 
deck,  a  Case  Control  deck,  and  a  Bulk  Data  deck.  Formats  for  all  cards 
used  here  are  defined  in  the  NASTRAN  User's  Manual  .  The  DMAP  ALTER  in¬ 
structions  contained  in  tite  Executive  Control  Deck  and  the  checkpoint 
dictionaries  apply  only  to  NASTRAN  Level  17. 

FLUID  MATRIX  CALCULATION 

To  obtain  a  series  of  fluid  H  and  L  matrices,  one  for  each  stage 
of  submersion,  one  static  stress  computer  run  is  made  for  each  stage. 

These  are  basically  static  stress  (Rigid  Format  1)  runs,  utilizing  CIUEXJ 
elements.  For  convenience,  all  grid  points  should  be  assigned  numbers 
higher  than  those  to  be  assigned  to  any  grid  points  of  the  propeller. 

The  element  and  grid  point  data  are  assembled  in  the  usual  way,  with  the 
following  exceptions  and  notations. 

1.  All  degrees  of  freedom  ul  all  grid  points  on  t  lie  outer  bound. irv, 
except  for  those  on  that  part  of  the  fluid-structure  interface  surface 
that  are  to  be  ii.  contact  with  the  propeller  during  llie  current  st.jge,  .ite 
constrained  wltli  SPC  cards. 

2.  All  but  one  translational  di-gree  ul  1  reedoin  ( c  oni|iotK'n  t  s  I,  ,  or 
3)  ot  all  other  grid  points  are  constrained. 


3.  Material  properties  are  input  on  a  MAT!  card.  The  value  ot 
the  shear  modulus  entered  is  the  reciprocal  of  the  mass  density  of 
water,  1/p,  for  the  present  case  in  which  water  contacts  only  one  side  of 
the  blade.  Young's  modulus  is  a  very  large  factor,  say  lU^*^,  times  the 
shear  modulus.  No  value  is  input  for  Poisson's  ratio,  so  that  the  NASTRAN 
program  will  compute  a  value. 

A.  Unit  pressure  is  applied  on  the  element  faces  jf  the  fluid  - 
structure  interface  which  are  to  be  in  contact  with  the  propeller  in 
the  present  stage.  Input  is  by  means  of  PL0AU3  cards.  The  unit  pressure 
can  be  either  positive  or  negative,  as  the  sign  will  not  affect  the  final 
result , 

5.  All  grid  points  in  contact  with  the  propeller  during  this  stage 
are  listed  on  ASETl  cards. 

The  Executive  Control  and  Case  Control  decks,  and  two  additional  cata¬ 
log  instructions  which  must  be  included  in  the  control  cards,  are  shown 
in  Eigure  A.l.  The  two  instructions  identify  two  permanent  files, 
UMATKIXSTEPI ,  and  LOAUSANUTABLESl ,  wliich  will  contain  the  H  [natrix  and 
other  information,  respectively. 

Eor  ttie  Second  and  succeeding  stages  ot  submersion,  the  only 
physical  change  is  that  a  new  set  of  grid  points  comes  into  contact  with 
the  propeller  blade.  Tills  set  will  generally  include  all  the  gridpoints 
included  in  the  previous  stage  of  submersion.  For  each  succeeding  run, 
then,  the  grid  points  newly  in  contact  with  the  propeller  must  be  added  to 
the  set  of  ASETl  points  and  removed  from  tlie  SP<,'  set  lor  which  all  degrees 
ot  t reedom  are  constrained.  The  element  faces  newly  in  lontact  with  the 
propeller  must  be  added  to  the  set  of  loaded  element  t -ices  ,  and  the  (^ATAl.iX. 
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cards  must  be  changed  to  assign  new  names  to  the  H  matrix  and  other 
data.  The  names  used  here  for  the  second  step  were  HMATK  1X.STEH2  and 
LUAUSANUTABLES2. 

The  printout  of  the  "PG"  matrix,  or  list  of  equivalent  forces,  is 
shown  in  Figure  lb.  It  was  necessary  to  obtain  this  matrix  before  con¬ 
straints  were  applied  in  order  to  obtain  all  three  components  tor  eac  li  ,  ■  'd 
point;  consequently,  tl>e  zeros  tor  all  other  unloaded  components  were 
also  printed  out.  Each  group  of  tliree  non-zero  numbers  represents  t he 
three  components  of  force  at  a  fluid  grid  point  in  contact  with  a  J■•ropeller 
grid  point.  The  numbers  are  given  in  tlie  order  in  whicli  the  fluid  crid 
points  are  numbered.  Tliese  values  are  to  be  entered  as  L*  ra.itrix  terms 
on  UMlG  cards  on  witich  ttie  matching  fluid  ;ind  propeller  interlace  grid 
points  will  also  be  listed.  The  procedure  tor  entering  these-  terms  on  I)Mlt. 
cards  is  best  understood  by  examining  tiie  example  to  be  presented  later. 

VIHRATION  ANAEYSIS 

The  Executive  Control  and  Case  Control  deiks  tor  the  vibration  aiialv- 
sis  of  a  submerged  propeller,  and  two  ATTACH  instructions  which  must  he 
included  with  the  control  cards,  are  g.iveii  in  ligiiri-  A.J.  Thi-  1  ile  naric-, 
on  tliese  two  eards  identity  the  d  m.itrix  .ind  .ithi-r  data  which  have  lain 
created  by  a  [irevious  run.  In  thi-,  example,  t  lie  ii.unei.  identifs  the  t  h  i  t -1  , 
or  1  ui  1  y-s  ubme  rged  ,  case.  The  AI.IEK  routine  in  the  l.xeiutice  lontii  1  Ui,. 
cont.iins  two  p.irtition  (I’AKI'N)  i  ns  t  nn- 1  1  oils  ,  one  ol  whiih  coul.iiic-  I  In 
matrix  naiiii- ,  I’VI.Ci.  That  name  was  arhilrarilv  hoseii  I-i  ideiilit-.  a 
pri  r  I  1  I  1  o  n  1  iig,  (  o  1  urn  11  matrix  o  t  the  same  name  t  -  -  hi  i  n  j  -  u  l  i  ii  tin-  hi  i  I  t  u  a  l  i  . 

The  hulk  data  is  set  Up  lias  i  ca  1  1  'c  as  a  p  1  opi  1  1  i  i  v  i  h  i  a  t  I  ■  -  n  lie;,  will 


.  1 


the  t  o  1  1  ow  1  u>;  a  dll  1  I  i  oils  : 


1.  All  flaid  points  in  contact  with  the  propeller  (in  this  fully- 
submerged  case,  all  fluid-structure  interface  points)  are  listed  as 
scalar  points  on  SPOINT  cards  in  this  propeller  run. 

2.  Constraints  on  the  zero  component  of  these  scalar  points  are 
input  on  SPC  cards . 

3.  Two  column  matrices,  PVECUi  and  PVEC3 ,  are  input.  In  tliis 
example  the  total  number  of  unconstrained  degrees  of  freedom  is  96x6, 

or  576  tor  the  propeller,  plus  40  for  the  scalar  points,  or  616.  Dimen¬ 
sions  of  tile  two  column  matrices,  therefore,  are  1x616.  All  entries  in 
PVECUI  are  zero,  and  all  except  the  last  4U  (assuming  the  SPOINT  numbers 
are  liigher  than  all  propeller  grid  point  numners)  of  the  PVEC3  entries  are 

Zero.  The  remaining  40  entries  are  1.0. 

T 

4.  The  L  matrix  terras  are  injiut  on  DMIC  cards.  In  this  case, 

40x3  terms  are  input. 

TRANSIENT  RESPONSE  ANALYSIS 

The  procedures  tor  preparing  the  additional  bulk  data  input  for  the 
vibration  analysis  are  essentially  the  same  as  those  for  the  transient 
response  case,  which  will  be  illustrated  in  more  detail. 

The  Executive  Control,  Case  Control,  and  Bulk  Data  decks,  and  the 
tour  additional  control  cards,  for  the  first  stage  of  submersion  of  the 
transient  response  case  are  shown  in  Figure  A. 4.  Of  the  four  control  cards, 
the  two  ATTACH  instructions  retrieve  the  aiiprojir i a t e  H  matrix  ior  this  ste|i, 
as  in  the  vibration  calculation.  The  REl^UEST  and  CATALOC  instructions 
ettect  the  storage  ot  the  checkpoint  till'  named  SfiHSTl’l,  required  for  the 
restart.  The  bulk  data  iniuit  consists  ot  the  same  data  as  tor  a  straight- 


torward  transient  response  analysis,  and  additional  SPOlN’l,  Sl’C,  DMIC,  and 


iM 


DMl  data  similar  to  that  for  the  vibration  in  water  case.  The  followinj^ 
comments  apply  to  the  additional  data; 

1.  All  fluid  grid  points  to  come  into  contact  with  the  propeller 
during  the  course  of  the  analysis  (not  just  those  in  contact  in  the  current 
stage  of  submersion)  must  be  listed  as  scalar  points  on  SHOINT  cards,  and 
also  on  SPC  cards  because  scalar  points  cannot  be  added  after  solution  of 
the  problem  has  begun. 

2.  PVtCDl  and  PVECl  matrices  have  the  same  dimensions  for  all  stages, 
in  tliis  example  lx  bib.  The  number  of  "l.U"  terms  in  PVECl,  however,  in 
this  case  18,  must  be  equal  to  the  number  of  fluid  points  in  contact  with 
the  propeller  for  the  current  stage,  and  the  placement  of  these  terms  with¬ 
in  the  "scalar  point  section"  of  the  PVECl  column  (tlie  last  41)  numbers) 
must  correspond  to  the  positions  of  the  18  contact  points  witliin  the 
"scalar  point  section."  All  degrees  of  freedom  of  propeller  grid  points 
and  scalar  points  are  represented  in  the  PVECl  column  in  numerical  order. 

In  this  example  tlie  fluid  points  were  numbered  so  that  the  portions  ol  t he 
blade  witli  the  highest  numbered  points  contact  tlie  water  first;  thus  <iny 

"i.O"  entries  in  PVECi,  or  PVEC2 ,  etc.,  are  the  last  terms  in  the  column. 

T 

3.  I'he  L  matrix  input  by  UMIC  cards  tor  tliis  stage  was  made  uji  t  roi.i 
the  corresponding  "PC"  matrix  printout  from  the  first  stage  iluid  run, 
given  in  Kigure  A. 2.  The  current  run  is  terminated  when  the  time  covered 
by  the  TSTEP  card,  that  is,  number  of  steps  times  lengtli  of  each  step  (in 
this  case  U. 91)93x10  ^  seconds),  has  elapsed. 

I'o  lonlinue  the  analysis  fur  the  second  stage  ut  submersion,  ,i  restart 
run  is  carried  out,  with  the  following  changes  to  the  deck: 


1.  The  two  ATTACH  instructions  and  the  CATALOG  instruction  described 
tor  the  first  run  are  replaced  with  the  following  control  cards: 

ATTACH, 0PTP,SUBSTP1,1D=CSRV. 

ATTACH , INPT , HMATR1XSTEP2 , 1 D=CSRV  . 

ATTACH, UTl,  LUADSANDTABLES2 , ID=CSRV  . 

CATALOG, NPTP, SUBSTP2 ,1D=CSRV. 

2.  The  checkpoint  dictionary  produced  by  the  previous  run  is  added  to 
the  Executive  Control  deck.  \lso,  PVECl,  whicli  appears  in  one  of  the 
PARTN  instructions,  is  ctianged  to  PVEt^d.  The  following  lines  are  added 
to  the  previous  ALTER  loutine  and  lelt  for  tlie  remaining  restart  runs: 

ALTER  UJ.lsJ 

t.HKPNT  kJPP,  K21Ji), >121)0,  B2DU,M1)L), KOI)  i> 

Al.TEK 

1.  Tile  rsi'El’  card,  TbrEP=2l,  in  tlie  Case  Control  deck,  is  replaced  with 
l'STEP  =  ..  . 

I.  Ila'  bu  i  K  liata  deck  consists  ot  t  lu'  lol  lowing  cards: 

I,  A  delete  (.  ,  )  card  which  eliminates  the  previous  DMIC  cards. 

T 

s.  y  1 1.  cards,  which  ri'preseut  tlie  1.  matrix  lor  stage  two. 

I  .  I'Ml  cards  lor  a  new  column  matrix,  PVEC2.  iJimensions  are  tlie 
sa::ie  as  tor  I’iECl,  but  the  number  ot  "i.O"  terms  is  changed  to 
represiMit  the  new  set  ot  points  in  contact  with  the  propeller. 

In  this  case,  t  In'  l.ist  2  b  tiTiiis  are  now  1.0. 
d.  A  iiew  iSii  r  bulk  liat.i  card,  identilied  as  TSTEP  number  22, 
to  correspond  to  the  l.ase  Control  iSTKP  card.  This  card 
should  cont.iin  the  lorrect  number  ot  time  stejis  so  that  tlu- 
total  time  delined  will  just  i'over  the  current  staga'  ol 

*  4 


submersion.  In  this  second  stage,  28  steps  were 
used . 

To  continue  with  tile  third  and  any  succeeding  steps,  the  procedures 
tor  the  second  stage  are  followed.  The  REQUEST  and  CATALOG  control  cards 
and  the  CHKRNT  Executive  Control  deck  card  are  not  required  in  the  final 
run . 


nast^an  executive  control  deck 


10  SES  PROPELLE'?  LUIO  lOEALlZATION 
APP  DISPLACEMENT 
SOL  1,0 

OIAG  1,6,1‘*,22 
TIME  60 
ALTER  97 

HATPRN  KAA, ,.,//$ 

OUTPUTl  ,,,,//C,N,-l  $ 

OUTPUTl  KAA,,,,//  $ 

JUMP  L3L7  S 

ALTER  ii2 

MATPRN  PG,ASET,,,//  5 

OUTPUTE  PG,EQEKIN, ASET, ,//C, N, 0/C, N.llS 

TA0PT  EQEXIN,.,.//  5 

JUMP  MEL  S 

ALTER  141 

label  MEL  B 

EQUIV  PG,PGG  S 

ALTER  143 

JUMP  L8L0EO  % 

ALTER  165 

JUMP  FINIS  « 

ENDALTER 

CENO 


CftSE  CONTROL  DECK 


title  3ES  PROPELLER  FLUID 
SUBTITLE  =  AODEO  MASS  TERMS 
SPC  =  11 
LOAD  =  1 

OISPLACcMENT  -  all 
makl Infs  =  5  ooo  oo 

BEGIN  BULK 

A.l  -  I’luid  kiiii  [’nrti.il  Fnpiit 


4  A 


S£S  FLJIO 

ftnOEO  “flSS  TEPHS 


TME  OFnSITy  of  TmIS  MfiTPIV  IS  5 . 7  PEPCENT 


ATTACH#INPT.HM*TRIXSTCP3.I0«CSSV. 
ATTACH.UTl ♦L0A0SAN0TABLrS3» 1D*CSRV. 


fcXFCUTlwt 


C  0  M  T  o  0 


DECK 


IJ  VI^sEATION  jta  P<OP  In 

APP  DISELACEM^nT 

SOL  3,C 

OIAG 

TIME  ol 

ALTEr  53 

OPO  dynamics , jPL , SIL ,USET/&PlO,SILD,USETD  EEO, EQOYN/ V,N, 

LUbET/W,N,LUSETCj/V,N,NOTFL/V,N,NOOLT/V,N,  NOPSDL/ V  ,  N,  NOFRL  /  V  . 
N,NONi.FT/W,N,NOTRL/»/  ,  N,  N0EE0/C,K,  /  V  ,  n,  nooe  $ 

SAVE  NOEEO.LuSETG  $ 

MTRXIN  CA5ECC,MA7PC0L,FQDYN,  » /K20PP  ,  M20P  P  ,  02  PP/ V  ,  N  ,  L  US<^T  0/ Y  » 

N,  N0K2  JPP/t/ I  N,N0M20PP/V  ,  N,  NO02PF  I 
PAPTN  K20PP, ,Ptf£C3/UNUSFD, ATMAT, ,/C,N,l  $ 

INPUTII  /, , , ,/C,N,-l  t 
INPUTII  /MMAT,,,,/  i 
SOLVE  HMAT , ATMA ( /HI AT/C,N, 1  i 
MPYAO  ATMAT,HIAT,MGG/MG&TEM/C,N,1  $ 

PARIN  '^C’iiTEH.PVECCl  ,/MGGNc  H»  ,  ,/C,N,-l/C  ,N,2/C,N  >6  $ 

PARAM  //C,N,MPY/V,N,POUM/C,N,l/C,H,-l  $ 

EQUIV  MGGNEWfMGG/ POUM  $ 

JJMP  oICK  S 

MATPKN  MGGNEKjMGG, , ,//  % 

LABEL  OICK  ^ 

ALTER  9E , 97 
ENOAlTE'-. 

CHKPNT  YtS 
CENT 


"-ASE  CONTROL  DECK 


TITLE=SEG  prop  VIBRATION 

SUBTITLE  =  PARTIALLY  SUBMERGED  PROPELLER 

LABEL  =  INVERSE  POKER  METHOD 

ECHO  =  both 

K2PP  =  STEF 

ME  Thr)n  =  4l 

SPC=1C 

OISPLACtMENT  =  All 
maxlines  =  SIODcJ 

BEGIN  BUL  K 


EiKure  A . i 


Propel  Ilt  Vibr.ition  Run  Partin]  Input 


Figure  A. 4  -  Propeller  Transient  Response  Run  Input  for  Stage  1 


REQUEST tNPTP.«PF. 
ATTACH,TNPTtHM*TRl*STEPl .tO»CSRV, 
ATTACH, UTl tLOAnSANOTAHLFSl tTO«CSRV. 
CATAL0G»NPTP,SUBSTP1 ,Tn«CSRV, 


NASTRan  EXECUTIVE  CONTROL  OFCX  ECHO 


TO  TRANSIENT  RESPONSE  OF  EES  PROPELLER  IN  WATER 

APP  niSPLACEMENT 

SOL  g,o 

OtAG 

TTME  60 

ALTER  72 

OPO  OYNAMTCS.hPL ,STL,USETTGPL0,5ILC,USET0,TFP00L,nLT, ,,NLFT,TRL , , 
EI}OYN/V,N,LUSET'/V,N,LUSfcTOy  V,N,N0TFL/V,K,N0nLT/V,N,N0PSOL/V, 
N,N0FRL/W ,N,NONLFT/V,N,NnTRL/ V,N,N0EEDFr,N,/V,N,N0UE  $ 

SAVE  LUSFirO,NODLT,NONLFT  .NOTRLjMOUE  t 

MTRYIN  rASEtC,H»TPOOL,eqDYN,,TPpoOLTK2OPP,M2Opp,02PP/V,N,LL'SFTO/'V,  N, 
NOt<?OPPXW  ,N,HO«aOPP/V,N,NOe2PP  t 
SAVE  NOXaepP.NOMenPp.NOpaPP  X 
PARTM  X?0PP, ,PVFC1/UNUSE0, ATNAT, ,/C,N,l  T 
INPUTTl  /,,  ,  ,/C,N  ,-l  S 
TNPUTTl  /HHat, ,, ,/  S 
SOLVE  HHAT, ATHATthIAT/C,N, I  X 
MPTAO  ATfi  AT,  VIAT,  mgC/NGOTEAFC,  N,  1  * 

PARTN  HSG’‘tH,PV<^Cei,/f«66NFW,,,/C,N,-l/C,N,2/C,N,6  !! 

EQOTV  liGGNEMtNGi-  * 

JtIHp  OTCK  * 

MATPRN  M&GNEH.MGr, ,,//  t 
LA5EL  OIC'<  ? 

ALTER 

ALTER  1?«,129 
ENOALTER 
CHKont  xfs 
CEND 
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Figure  A. 4  (Continued) 
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